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Abstract 

Evidence for the onset of deconfinement in central Pb+Pb collisions was reported by NA49 in fixed target measure¬ 
ments at beam momentum 30A GeV/c. This observation motivated the NA61/SHINE program started in 2009 at the 
CERN SPS, which, in particular, aims to study properties of the onset of deconfinement by measurements of hadron 
production in proton-proton, proton-nucleus and nucleus-nucleus collisions. 

This contribution presents spectra of charged pions produced in p-rp interactions and ^Be-r^Be collisions at 20A- 
158A GeV/c beam momentum. The NA61/SHINE results are compared with the corresponding NA49 data from 
central Pb-hPb collisions at the same beam momenta per nucleon. 
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1. Introduction 

The NA61/SHINE experiment at the CERN SPS con¬ 
ducts a broad programme of hadron production mea¬ 
surement in proton-proton, nucleus-nucleus, proton- 
nucleus and pion-nucleus collisions. This paper focuses 
on measurements of collisions of protons and nuclei at 
different projectile momenta. A two-dimensional scan 
is performed in beam energy from 13A to 158A GeV/c 
and in system size from p-hp up to Xe-rLa. It is planned 
to repeat and extend measurements of the NA49 ex¬ 
periment on Pb-l-Pb collsions (nia with the upgraded 
NA61/SHINE setup. Eigure[2 visualises the scan. 

The scan aims to search for the critical point of the 
quark-gluon plasma phase transition and study the on¬ 
set of deconfinement in collisions of light nuclei. The 
results presented in this paper relate to the latter topic. 

This paper presents double differential spectra of neg¬ 
atively charged pions produced in p-rp interactions at 
20, 31, 40, 80 and 158 GeV/c (published O) and in 
^Be-r^Be interactions at 40A, 75A and 150A GeV/c (pre¬ 
liminary), and positively charged pion spectra produced 


in p-hp interactions at 40, 80 and 158 GeV/c (prelimi¬ 
nary). 


2. Experimental set-up 

2.7. Beams and targets 

The experiment is located in the H2 beamline of the 
CERN North area. Secondary protons were selected 
with the desired beam momentum and directed on a 
20 cm-long cell filled with liquid hydrogen served as 
the proton target. 

The beryllium beam was produced by fragmentation 
of the primary Pb beam on a beryllium target. The H2 
beam line was operated as a fragment separator. In or¬ 
der to ensure purity of the beam long-living ^Be nuclei 
were selected, as neighbouring isotopes: ^Be and ^Be 
are unstable and ^^Be cannot be separated sufficiently 
from ^Be. The interaction target was made of stable 
^Be. 
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method allows to derive the 7i~ spectra in a broad phase- 
space region. The dE/dx identification method is does 
not work at /? < 4 GeV/c where the Bethe-Bloch curves 
cross and at high momenta, where the statistics is insuf¬ 
ficient to fit the data. 

The h“ method was used to derive the 7r~ spec¬ 
tra in the full geometrical acceptance in p-hp interac¬ 
tions at 20-158 GeV/c and in Be-hBe collisions at 40- 
150A GeV/c. The dE/dx method was used to derive the 
7r“ and spectra in p-rp interactions at 40-158 GeV/c 
in a more region. The results were corrected for the de¬ 
tector acceptance and efficiency and other experimental 
effects. 

3.2. Centrality selection with PSD 


Figure 1: NA61/SHINE program of ion-ion collisions: two- 

dimensional scan in the system size (from p-i-p to Pb-i-Pb) and the 
projectile momentum (from 13A to 158A GeV/c). Each small box cor¬ 
responds to 2-5 million of events; large boxes are 50 million. The 
green boxes show the data collected as of summer 2014, and the blue 
shade marks the data analysed in this paper. 

2.2. NA61ISHINE detector 

Figure shows a schematic of the NA61/SHINE de¬ 
tector (upgrade of NA49 apparatus). Well-aligned beam 
particles are detected in a set of scintillation counters 
(SI, S2, VO, VI, Vr). An additional counter S4, used 
in anti-coincidence, was placed behind the target on the 
extrapolated beam trajectory. The beam is identified 
by two Cherenkov detectors: CEDAR and Threshold 
Cherenkov Counter (THC). The Be beam is also iden¬ 
tified by the A detector (time of flight) and the Z de¬ 
tector (Cherenkov). The Projectile Spectator Detector 
(PSD) calorimeter on the extrapolated Be beam trajec¬ 
tory measures the non-interacting nucleons (spectators) 
of the projectile nuclei. Three Beam Position Detectors 
(BPD) measure the transverse coordinates of the inter¬ 
action point. Produced charged particles are tracked in 
five Time Projection Chambers: two VTPCs placed in¬ 
side magnets, two MTPCs and the GAP TPC. Time of 
Flight (ToF) walls are placed behind the MTPCs. 


Centrality of the ^Be-f^Be collisions was determined 
based on the forward energy measured with the PSD. 
Four centrality regions were defined: 0-5-10-15-25%. 

Two effects shift the rapidity spectra in opposite di¬ 
rections: 

• mass asymmetry between the projectile and the tar¬ 
get biases the spectra towards negative rapidity, 

• centrality selection with the PSD biases the spectra 
towards positive rapidity. 

These effects cancel each other partially, however their 
impact on the results is still under investigation. For this 
reason only the mid-rapidity spectra will be discussed 
for Be-fBe collisions. 


4. Results 

4.1. Double differential spectra 

Double differential tt spectra are shown in Figs. 
and|^ The 7r~ spectra were derived with the h“ method 
in a large region of phase-space. The phase-space avail¬ 
able for the TT^ spectra obtained with the dE jdx method 
is the largest at 158 GeV/c. The missing region at low y 
and low /7 t will be covered by the ToF analysis. 


3. Analysis method 

3.1. Derivation of spectra 

The majority of the produced negatively charged 
hadrons are pions. A small (<10%) contribution of 
other hadrons, mostly p and K“ can be reliably removed 
based on model predictions. This so-called h“ analysis 


4.2. Comparison with the available data 

Figure shows integrated rapidity spectra of 7i~ and 
TT^ mesons in p-fp interactions. Within the systematic 
uncertainty the results of the h“ and the dE jdx methods 
agree. They also agree with the measurements of other 
experiments. 
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Figure 2: Schematic of the NA61/SHINE detector (horizontal cut, not to scale). The beam line set-up is magnified at the bottom of the figure. 
The spectrometer is formed by five Time Projection Chambers (TPC, blue). VTPCs are located inside VTX magnets. Behind the MTPCs there are 
three Time Of Flight (ToF) planes and the Projectile Spectator Detector (PSD) calorimeter. 


4.3. 71 transverse mass spectra at mid-rapidity 

The transverse mass (mx = ^Pj + m^) spectra of 
7i~ in p+p interactions and Be+Be collisions at 40A- 
150A GeV/c are shown in ¥ig.^(top). The spectra are 
fitted with an exponential function 

1 d^ ( mT\ ... 

mx dmx V T J 

in 0.2 < mx < 0.7 GeV/c^. The fitted inverse slope 
parameter T is plotted against the number of wounded 
nucleons Nw (number of nucleons interacting inelasti- 
cally) in Vig.^ihottom). T is seen to increases mono- 
tonically with Nw 

Figure]^ shows the energy dependence of the inverse 
slope parameter T. For Pb+Pb the T parameter is larger 
by 10-20 MeV than for p+p interactions. For Be+Be at 
40A GeV/c T is the same as for p+p. At higher energies 
T increases by more than 10 MeV for Be+Be, while it is 
almost constant for p+p. This suggests the onset of col¬ 
lective effects for Be+Be collisions above 40A GeV/c. 

Figure shows the ratio of the transverse mass spec¬ 
tra of 7i~ mesons produced at mid-rapidity in central 
Pb+Pb collisions and inelastic p+p interactions at the 
same collision energy per nucleon. The spectra were 
normalised to unity before dividing. The ratio ex¬ 
ceeds unity for mx - m;^ < 0-1 GeV/c^ and mx - m;^ > 


0.5 GeV/c^, while it is below unity in 0.1 < mx - m;^ < 
0.5 GeV/c^, demonstrating different shapes of the spec¬ 
tra. The ratio shows weak dependence on the collision 
energy. 

Comparison of p+p interactions and collisions of Pb 
nuclei, composed of protons and neutrons, requires tak¬ 
ing into account the isospin symmetry. Full interpreta¬ 
tion of the results requires spectra in p+p derived in 
a large region of phase-space. 


4.4. tT rapidity spectra 

Figure shows the integrated rapidity spectra of 7i~. 
The spectra are described well by the sum of two sym¬ 
metrically displaced Gaussian distributions: 
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Figure [T^ shows the ratio of the n~ rapidity spectra 
normalised to unity in central Pb+Pb collisions and in¬ 
elastic p+p interactions. Rapidity was normalised to 
beam rapidity ybeam- The distributions in Pb+Pb are 
broader as at y > 0.7 the ratios exceed unity by 10- 
20%. The energy dependence is weak. Spectra of 
in full phase-space are needed to understand the differ¬ 
ences. 
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Figure 3: Double differential spectra of n in p+p interactions at 20, 31, 40, 80 and 158 GeV/c (first and second row) and spectra in p+p 
interactions at 40, 80 and 158 GeV/c (second row). 
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Figure 4: Double differential spectra of n in ^Be-i-^Be collisions at 40A, 75A and 150A GeV/c in four centrality classes. 
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Figure 5: The NA61/SHINE integrated rapidity spectra of n~ and mesons in inelastic p+p interactions at 31 and 158 GeV/c obtained with the 
h“ (left, middle) and dE/dx methods (middle, right), compared with results from MIRABELLE at 32 GeV/c (4) and NA49 at 158 GeV/c (5). The 
shaded bands show the systematic uncertainty of the NA61/SHINE data. 
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Figure 6: Top: n transverse mass spectra at mid-rapidity (0 < y < 0.2) for inelastic p-i-p interactions, Be-i-Be collisions in four centrality regions 
and central Pb-i-Pb collisions [T] at 40A (left), 75/80A (middle) and 150/158A GeV/c (right). The spectra are fitted with an exponential function 
Eq.Q in the regions of 0.2 < mx < 0.7 GeV/c^, shown by thick solid lines. Bottom: fitted values of the inverse slope parameter T as a function of 
the number of wounded nucleons Nw for inelastic p-i-p interactions and Be-i-Be and Pb-i-Pb collisions. 
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Figure 7: Inverse slope parameter T of the n transverse mass dis¬ 
tributions in inelastic p-i-p and central Be-i-Be and Pb-i-Pb (QIZ) col¬ 
lisions. The T parameter was fitted in the region of 0.2 < mx < 
0.7 GeV/c^ 


Figure 9: Integrated rapidity spectra of n~ produced in inelastic p-i-p 
interactions at 20-158 GeV/c. The open points are reflected with re¬ 
spect to y = 0. The shaded bands show the systematic uncertainty. 
The curves show function Eq. ^ fitted to the data. 
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Figure 8: Ratio of the normalised n transverse mass spectra at mid¬ 
rapidity produced in central Pb-i-Pb collisions EM and inelastic p-i-p 
interactions at the same collision energy per nucleon. The coloured 
bands represent the systematic uncertainty of the p-i-p data. 


Figure 10: Ratio of normalised tt integrated rapidity spectra of tt in 
central Pb-i-Pb collisions dill] and inelastic p-i-p interactions at 20A- 
158A GeV/c. The shaded bands show the systematic uncertainty of the 
p-i-p data. 
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Figure 11: Mean multiplicity of all pions per wounded nucleon 
produced in inelastic p+p interactions and central Pb+Pb mm 
(Au+Au muill]) collisions versus Fermi energy F ^ s^^^^ (9j- The 
vertical lines show the total uncertainty. The open points show results 
from the other experiments, mostly using bubble chambers (5][Tol. 

4.5. Total pion multiplicity 

The total 7i~ multiplicity was calculated by inte¬ 
grating the rapidity spectra. Known phenomenolog¬ 
ical relations between multiplicities of n of different 
charges ifTTl allowed to calculate the multiplicity of n 
of all charges. Figure pT] shows the energy dependence 
of the pion multiplicity (tt) divided by the number of 
wounded nucleons in inelastic p-fp and central Pb-fPb 
collisions. The NA61/SHINE p-fp results agree with re¬ 
sults from other experiments. The multiplicity increases 
linearly with the Fermi energy F for p-fp interactions. 
For central Pb-fPb collisions the increase is faster in the 
SPS energy region. This supported the conclusion that 
the onset of deconfinement in central Pb-fPb collisions 
occurs near 30A GeV/c lO. The NA61/SHINE results 
demonstrate the capability of the experiment to study 
the onset of deconfinement in light ion collisions. 

5. Conclusions 

This paper presents the n~ spectra in full phase-space 
in p-fp interactions at 20-158 GeV/c and in ^Be-f^Be 
collisions at 40A-150A GeV/c, as well as spec¬ 
tra in limited phase-space in p-fp interactions at 40- 
158 GeV/c. The results contribute to the NA61/SHINE 
programme of study of onset of deconfinement. 

The inverse slope of the 7:~ transverse mass distribu¬ 
tion in Be-fBe collisions increases with energy above 


40A GeV/c, while it is almost constant in p-fp. This sug¬ 
gests the onset of collective effects in Be-fBe collisions 
above 40A GeV/c. Spectra of 7:~ in p-fp interactions and 
in central Pb-hPb collisions differ significantly. 

Full interpretation of the presented results requires 
further investigation of the forward-backward asymme¬ 
try in Be-rBe collisions. Spectra of in full phase- 
space are required to understand the role of isospin ef¬ 
fects. 
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